arXiv:1502.01724vl [nucl-th] 5 Feb 2015 


PROCEEDINGS 

SCIENCE 



Direct photon production at RHiC and LHC energies 


O. Linnyk* 

Institutfiir Theoretische Physik, Universitat Giefien, 35392 Giefien, Germany 

E-mail: |linnyk@f ias . uni-f rankf urt. de 

E. L. Bratkovskaya 

Institutfiir Theoretische Physik, Universitat Frankfurt, 60438 Frankfurt am Main, Germany 

W. Gassing 

Institutfiir Theoretische Physik, Universitat Giefien, 35392 Giefien, Germany 

Direct photon spectra and elliptic flow V 2 in heavy-ion collisions at RHIC and LHC energies 
are investigated within a relativistic transport approach incorporating both hadronic and partonic 
phases - the Parton-Hadron-String Dynamics (PHSD). The results suggest that a large V 2 of the 
direct photons - as observed by the PHENIX Collaboration - signals a significant contribution 
of photons produced in interactions of secondary mesons and baryons in the late stages of the 
collision. In order to further differentiate the origin of the direct photon azimuthal asymmetry, we 
compare our predictions for the centrality dependence of the direct photon spectra to the recent 
measurements by the PHENIX Collaboration and provide predictions for Pb + Pb collisions at 
LHC energies with respect to the direct photon spectra and V 2 (pr) for 0-40% centrality. 


XXII International Baldin Seminar on High Energy Physics Problems 
15-20 September, 2014 
JINR, Dubna, Russia 


* Speaker. 


(c) Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 


http://pos.sissa.it/ 











Direct photons at RHIC and LHC 


O. Linnyk 


1. Introduction 


Real and virtual photons are powerful tools to probe matter under extreme conditions as cre¬ 
ated in heavy-ion collisions at relativistic energies, since electromagnetic radiation is emitted over 
the whole collision evolution. The first several fm/c of the collision are particularly interesting be¬ 
cause of the high energy densities reached (well above the expected phase transition region to the 
deconfined phase). The photons interact only electromagnetically and thus escape to the detector 
undistorted through the dense and strongly-interacting initially produced medium. Their differen¬ 
tial spectra and elliptic flow carry the information on the properties of the matter produced to the 
detector. 

On the other hand, the measured photons provide a time-integrated picture of the heavy-ion 
collision dynamics and are emitted from every moving charge - partons or hadrons. Therefore, a 
multitude of photon sources has to be differentiated in order to access the signal of interest. The 
dominant contributions to the inclusive photon production are the decays of mesons, mainly pions, 
eta- and omega-mesons. The PHENIX and ALICE Collaborations subtract the “decay photons" 
from the inclusive photon spectrum using a cocktail calculation ||i|, ^ and obtain the “direct" pho¬ 
tons. 


In particular the direct photons at low transverse momentum {pj < 3 GeV) are expected to be 
dominated by the "thermal" sources, i.e. the radiation from the strongly interacting Quark-Gluon- 
Plasma (sQGP) llll] as well as secondary meson-i-meson and meson-i-baryon interactions [§, 
These partonic and hadronic channels have been studied within PHSD in detail in Refs. [|0il at 
Relativistic-Heavy-Ion-Collider (RHIC) energies. It was found that the partonic channels constitute 
up to half of the observed direct photon spectrum for very central collisions. Other theoretical 
calculations find a dominant contribution of the photons produced in the QGP to the direct photon 
spectrum 

The low-pT direct photons probe not only the temperature |[I], ^ |T^ of the produced QCD- 
matter, but also its (transport) properties, for instance, the sheer viscosity. Using the direct photon 
elliptic flow V 2 (a measure of the azimuthal asymmetry in the photon distribution) as a viscometer 
was first suggested by Dusling et al. in Ref. [13]; this idea was later supported by the calculations in 
Refs, [ll], 14]. It was also suggested that the photon spectra and V 2 are sensitive to the collective 

directed flow of the system [15, and to the asymmetry induced by the strong magnetic field 
(flash) in the very early stage of the collision [|^, 

However, the recent observation by the PHENIX Collaboration [|T]] that the elliptic flow V 2 {pt) 
of ’direct photons’ produced in minimal bias Au-i-Au collisions at y/s^N = 200 GeV is comparable 
to that of the produced pions was a surprise and in contrast to the theoretical expectations and 
predictions. Indeed, the photons produced by partonic interactions in the quark-gluon plasma phase 
have not been expected to show considerable flow because they are dominated by the emission in 
the initial phase before the elliptic flow fully develops. 

In Ref. [^ we have applied the PHSD approach to photon production in Au-i-Au collisions 
at y/SNN = 200 GeV and studied the transverse momentum spectrum and the elliptic flow V 2 of 
photons from hadronic and partonic production channels. A microscopic description of the full 
collision evolution is done by the covariant off-shell transport PHSD. The degrees of freedom in 
the partonic and hadronic phases are strongly interacting dynamical quasi-particles and off-shell 
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hadrons, respectively. Further, it was found in Ref.0 that the PHSD calculations reproduce the 
transverse momentum spectrum of direct photons as measured by the PHENIX Collaboration in 


Refs. [ ]19i |20t ]. The centrality dependence of the thermal photon yield in PHSD was predicted 
to be ~ the exponent a = 1.5, which is in a good agreement with the most recent 

measurement of a = 1.48 ±0.08 ± 0.04 by the PHENIX Collaboration [^T|]. Eurthermore, the 
PHSD also described the data on the elliptic flow of inclusive and direct photons. The strong V 2 
of direct photons - which is comparable to the hadronic V 2 - in PHSD is attributed to hadronic 
channels, i.e. to meson binary reactions which are not subtracted in the data. As sources for 
photon production, we incorporated the interactions of off-shell quarks and gluons in the strongly 
interacting quark-gluon plasma (sQGP) {q + q ^ g + J and q{q) ± g —>• q{q) + j), the decays of 
hadrons (tt —> 7± 7, Tj — ^ 7± 7, ft) —)■ tt ± 7, —)• p ± 7, 0 — ^ ± 7, ai —^ tt ± 7) as well as their 

interactions (7r±7r—)-p±7, p±7r—)'7r±7, meson-meson bremsstrahlung m + m ^ m + m + y), 
meson-baryon bremsstrahlung (m±R —)■ m±R± 7) and the two-to-two meson-i-baryon interactions 
(p ±p —> 7± p/n and p ±n —^ 7± p/n) . 

The photon production via bremsstrahlung in meson-meson and meson-baryon elastic colli¬ 
sions was found to be an important source for the direct photon spectra and elliptic flow simultane¬ 
ously lH, 0], where for the calculation of the photon bremsstrahlung from all elastic meson-meson 
and meson-baryon scatterings m\ ± m 2 , which occur during the heavy-ion collisions (including 
nij = n,ri,K,K,K^,K*,K*,K*^,7]', 0 i),p,(j),ai),'we have been applying the soft photon approxima¬ 
tion. Therefore the resulting yield of the bremsstrahlung photons depended on the model assump¬ 
tions such as (i) the cross section for the meson-meson elastic scattering (we assumed 10 mb for 
all meson species), (ii) incoherence of the individual scatterings and (iii) the soft photon approxi¬ 
mation (i.e. low photon energy and low ^/s of the collision). The adequacy of the SPA assumption 


has been studied in Ref. |22] and a theoretical uncertainty of up to a factor of 2 was found. 

The results of our calculations so far have been compared to the data from RHIC. Additionally, 
here we will provide calculations for the photon production in Pb ± Pb collisions as the energy 
of y/sNM = 2.76 TeV. Since the preliminary data of the AEICE Collaboration indicate a 

significant direct photon signal at low pj with a large elliptic flow at EHC energies, a differential 
comparison of our calculations with the final data will be mandatory. 

2. PHSD 

To address the photon production in a hot and dense medium - as created in heavy-ion colli¬ 
sions - we employ an up-to-date relativistic transport model, i.e. the Parton Hadron String Dynam¬ 


ics [ |24{ ] (PHSD) that incorporates the explicit partonic phase in the early reaction phase. Within 
PHSD, one solves generalized transport equations on the basis of the off-shell Kadanoff-Baym 
equations for Greens functions in phase-space representation (in first order gradient expansion, be¬ 
yond the quasiparticle approximation). The approach consistently describes the full evolution of a 
relativistic heavy-ion collision from the initial hard scatterings and string formation through the dy¬ 
namical deconfinement phase transition to the quark-gluon plasma (QGP) as well as hadronization 
and to the subsequent interactions in the hadronic phase. In the hadronic sector PHSD is equivalent 
to the Hadron-String-Dynamics (HSD) transport approach [ ^ ] that has been used for the descrip¬ 
tion of pA and AA collisions from SIS to RHIC energies and has lead to a fair reproduction of 
hadron abundances, rapidity distributions and transverse momentum spectra. The description of 
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quarks and gluons in PHSD is based on a dynamical quasiparticle model for partons matched to 
reproduce lattice QCD results in thermodynamic equilibrium (DQPM). The DQPM describes QCD 
properties in terms of single-particle Green’s functions (in the sense of a two-particle irreducible 
approach) and leads to the notion of the constituents of the sQGP being effective quasiparticles, 
which are massive and have broad spectral functions (due to large interaction rates). The transi¬ 
tion from partonic to hadronic degrees of freedom in PHSD is described by covariant transition 
rates for the fusion of quark-antiquark pairs to mesonic resonances or three quarks (antiquarks) 
to baryonic states. The PHSD transport approach provides a good description of the data on bulk 
properties as well as hard probes |27] for a wide range of energies up to the LHC. 


3. Direct photon production 

We consider the following sources of direct photons: 

1) Photons radiated by quarks in the interaction with other quarks and gluons in the two-to-two 
reactions: 


q+q^g+y, 

q/q+g^ q/q+y- 

The implementation of the photon production by the quark and gluon interactions in the PHSD 
is based on the off-shell cross sections for the interaction of massive dynamical quasi-particles as 
described in [^, In addition, photon production in the bremsstrahlung reactions q + q/g ^ 
q + q/g + y should be incorporated. 

2) All colliding hadronic charges (meson, baryons) can also radiate photons by the bremsstrahlung 
process: 


m + m^m + m + Y^ 
m + B^m + B + y. 


(3.1) 

(3.2) 


The processes (^H]) have been calculated within the PHSD in Refs. [^, |^, while the m + B 
bremsstrahlung ( |3.2| ) reactions have been added in Ref. 0. The implementation of photon bremsstrahlung 
from hadronic reactions in transport approaches so far has been based on the ’soft photon’ approx¬ 


imation (SPA) []30|], which relies on the assumption that the radiation from internal lines is negli¬ 
gible and the strong interaction vertex is on-shell; this is valid only at low energy (and pj) of the 
produced photon. 

3) Additionally, the photons can be produced in specific binary hadronic collisions. We con¬ 
sider the direct photon production in the following 2—^2 meson-i-meson collisions 


n + n^ p + y, 
n + p ^n + y, 


accounting for all possible charge combinations. The implementation of these reactions has been 
described in Refs. |^, ^]. 
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Figure 1: PHSD results for the spectrum of direct 
photons produced in 0-40% most central Auh-Au col¬ 
lisions at y/s^ = 200 GeV as a function of the trans¬ 
verse momentum pr at mid-rapidity |y| < 0.5. The 
data of the PHENIX Collaboration are taken from 
Ref. @. 



Figure 2: Elliptic flow V 2 versus transverse momen¬ 
tum pt for the direct photons produced in the mini¬ 
mal bias Auh-Au collisions at ^/s^ = 200 GeV calcu¬ 
lated within the PHSD (solid red line); the blue band 
reflects the uncertainty in the modeling of the cross 
sections for the individual channels. The data of the 
PHENIX Collaboration are from Ref. [|^. 


4. Results 

The results for the direct photon spectrum as a sum of partonic as well as hadronic sources 
for the photons produced in 0-40% central Au-i-Au collisions at y/sm = 200 GeV is presented 
in Fig. H as a function of the transverse momentum px at mid-rapidity |y| < 0.5. The calculated 
channel decomposition of the spectrum is presented in Fig. |I| by the lines of various styles. The 
measured transverse momentum spectrum dNjdpj (given by the filled circles) is reproduced well 
by the sum of partonic and hadronic sources (red solid line). 

We find that the radiation from the sQGP constitutes slightly less than half of the observed 
number of photons. The radiation from hadrons and their interaction - which are not measured 
separately so far - give a considerable contribution at low transverse momentum. The dominant 
hadronic sources are the meson decays and the meson-meson bremsstrahlung. While the former 
(e.g. the decays of ft), T]’, 0 and a\ mesons) can be subtracted from the spectra once the mesonic 
yields are determined independently by experiment, the reactions n + p^n + y, n + ’K^p + y, 
p+p/n^n/p + y and the meson-meson and meson-baryon bremsstrahlung can be separated from 
the partonic sources only using theoretical models. 

The azimuthal momentum distribution of the emitted particles is commonly expressed in the 
form of Fourier series as 


d'^N 


d^N 


InprdpTdy 


1 + Y,2v„{pT)cOs[n{\l/-'i'n)] 


n—l 


(4.1) 


where v„ is the magnitude of the nth order harmonic term relative to the angle of the initial- 
state spatial plane of symmetry 'F„. One should take into account event-by-event fluctuations 
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|GeV/c| |GeV/c| 

Figure 3: Contribution of the photon production in the two-to-two p+nucleon interaction (orange dash 
lines) to the total direct photon spectra (red lines) at the top RHIC energy of ^Jsnn = 200 GeV at differ¬ 
ent centralities. The dominant sources are the photons from the QGP and from the hadronic two-to-three 
bremsstrahlung processes. The theory lines - except the p-rnucleon channel - are taken from the PHSD 
predictions published in . The PHENIX data are from Ref. [ ^0| , pl| ] . 


with respect to the event plane We calculate the V 3 coefficients with respect to *^3 as 

V 3 {*P 3 } = (cos(3[t/r —'P 3 ]))/Res(*P 3 ). The event plane angle *^3 and its resolution Res(*P 3 ) are 
calculated as described in Ref. I0- We recall that the second flow coefficienl V 2 carries informa¬ 
tion on fhe inferacfion sfrength - and fhus on fhe sfafe of mailer and ils properlies - al Ihe space-lime 
poinl, from which fhe measured particles are emilled. 

Aboul a decade ago, Ihe WA98 Collaboration has measured Ihe elliptic flow V 2 of photons 
produced in Pb+Pb collisions al Ihe beam energy of Ebeam = 158 AGeV [|^], and il was found lhal 
Ihe V 2 ( 7 '”'^^) of Ihe low-transverse-momenlum inclusive photons was equal to Ihe V 2 ( 7 ^) of pions 
wilhin Ihe experimenlal uncerlainlies. This observation lead to Ihe conclusion lhal eilher (Scenario 
a:) Ihe conlribulion of Ihe direcl photons to Ihe inclusive ones is negligible in comparison to Ihe 
decay photons, mainly Ihe 71 ^ decay producls, or (Scenario 2:) Ihe elliptic flow of Ihe direcl photons 
is comparable in magnilude to Ihe V 2 ( 7 ""^ 0 > V 2 ( 7 ^'^“^) and V 2 ( 7 r). 

However, in view of Ihe WA98 measuremenl of Ihe direcl photon speclrum, which we de¬ 
scribed above, Ihere is a significanl finite yield of direcl photons al low Iransverse momenlum. 
Thus Ihe scenario 1 can be ruled out Furthermore, Ihe observed direcl photons of low pj musl 
have a significanl elliptic anisolropy V 2 of Ihe same order of magnilude as Ihe hadronic flow. This 
scenario poinls towards hadronic sources of Ihe direcl photons. Thus, Ihe inlerprelalion [ p^ , ^ of 
Ihe low-pr direcl photon yield measured by WA98 - as dominanlly produced by Ihe bremsslrahlung 
process in Ihe mesonic collisions TT + TI—)-7r + 7r + 7 -isin accord also wilh Ihe dala on Ihe photon 
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Figure 4: Yield of direct photons in Pb+Pb collisions at the invariant energy = 2.76 TeV for 0-40% 
centrality as predicted within the PHSD. 


elliptic flow V 2 ( 7 ""^^). 

Let us note that the same conclusions apply also to the most recent studies of the photon elliptic 
flow at RHIC and LHC. The PHENIX and ALICE Collaborations have measured the inclusive 
photon V 2 and found that at low transverse momenta it is comparable to the V 2 {pt) of decay photons 
as calculated in cocktail simulations based on the known mesonic V 2 {pt)- Therefore, (a) either the 
yield of the direct photons to the inclusive ones is not statistically significant in comparison to 
the decay photons or (b) the elliptic flow of the direct photons must be as large as V 2 ( 7 ^'^“'’) and 

In the PHSD, we calculate the direct photon V 2 ( 7 ^"^) by building the weighted sum of the 
channels, which are not subtracted by the data-driven methods, as follows: the photons from the 
quark-gluon plasma, from the initial hard parton collisions (pQCD photons), from the decays of 
short-living resonances (a-meson, 0-meson, A-baryon), from the 2 —)• 2 channels (tt + p —)• tt + 7 , 
71 + 71 —)• p + 7 ), and from the bremsstrahlung in the elastic meson-i-meson and meson-i-baryon 
collisions {m + m^m + m + Y,m-\-B^m + B + Y). We calculate the direct photon V 2 (in PHSD) 
by summing up the elliptic flow of the individual channels contributing to the direct photons, using 
their contributions to the spectrum as the relative pT'-dependent weights, w,(p 7 ’), i.e. 


V2(/"') = L^2(7'VKP7’) 
i 


LiV2{Y)Ni{pT) 

ZMpt) 


(4.2) 


The results for the elliptic flow V 2 {pt) of direct photons produced in Au +Au collisions at the top 
RHIC energy are shown in Pig. ^ According to our calculations of the direct photon spectra, almost 
a half of the direct photons measured by PHENIX stems from the collisions of quarks and gluons 
in the deconfined medium created in the initial phase of the collision. The photons produced in the 
QGP carry a very small V 2 and lead to an overall direct photon V 2 about a factor of 2 below the pion 
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Figure 5: Left: Elliptic flow V 2 versus transverse momentum pj for the inclusive photons produced in 0- 
40% central Pb+Pb collisions at = 2.76 TeV as predicted by the PHSD (solid red line); the blue error 
band reflects the bnite statistics and the uncertainty in the modeling of the cross sections for the individual 
channels. Right: Elliptic flow vi versus transverse momentum pr for the direct photons produced in 0-40% 
central Pb-tPh collisions at ^/snn = 2.76 TeV as predicted by the PHSD (solid red line); the blue error band 
is dominated by the uncertainty in the modeling of the cross sections for the individual channels. 


V2{7 i) even though the other channels in the sum ( |4.2[ ) have large elliptic flow coefficients V 2 of the 
order of V2(7r) (cf. Ref. Ip). 

Indeed, the parton collisions - producing photons in the QGP - take place throughout the 
evolution of the collision but the collision rate falls rapidly with time and thus the production of 
photons from the QGP is dominated by the early times (cf. Fig. 7 in Ref. p). As a consequence, the 
elliptic flow ‘picked up’ by the photons from the parent parton collisions saturates after about 5 fm/c 
and reaches a relatively low value of about 0.02, only. We note that a delayed production of charges 
from the strong gluon fields (‘glasma’ l^]) might shift the QGP photon production to somewhat 
later times when the elliptic flow is built up more. However, we cannot quantitatively answer 
whether the additional evolution in the pre-plasma state could generate considerable additional V 2 . 

In addition, the centrality dependence of the direct photon spectra and flow has been investi¬ 
gated. The recent measurements by the PHENIX Collaboration [ H | confirm the predictions within 
the PHSD from Ref. [Q]. We present our predictions and the data in Fig. |3[ The centrality depen¬ 
dence of the integrated thermal photon yield in PHSD was found to scale as with the exponent 
a= 1.5, which is in a good agreement with the most recent measurement of a = 1.48 ± 0.08 ± 0.04 
by the PHENIX Collaboration [H]. 

Centrality dependence of the direct photon elliptic flow V 2 was also calculated within the PHSD 
in Ref. 0 and confirmed (within the error bars) by the PHENIX Collaboration measurements in 
Ref. [ ^0| , ^T| ]. Thus the observed centrality dependence of the elliptic flow is in agreement with the 
interpretation of the direct photons having a hadronic origin (in particular from the bremsstrahlung 
in meson-i-meson and meson-i-baryon collisions), which is stronger in more peripheral collisions. 

In Fig. ^ we show the calculated direct photon yield in Pb-i-Pb collisions at the invariant energy 
Hsmn = 2.76 TeV for 0-40% centrality. Comparing the theoretical predictions to the preliminary 
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data of the ALICE Collaboration from Ref. 0, we find an overall agreement with the data within 
about a factor of 2 in the range of transverse momenta pj from 1 to 4 GeV. On the other hand, the 
calculations tend to underestimate the preliminary data in the losN-pj region and to underestimate 
slightly the highest-pr points. One cannot exclude the existence of another yet unknown source of 
direct photons at low pj. However, the significance of fhe comparison is nof robusf unfill fhe final 
dafa will be available. 

We finally presenf our predictions for fhe elliptic flow of inclusive and direcf phofons produced 
in Pb + Pb collisions af fhe energy of y/¥m = 2.76 TeV af fhe LHC wifhin fhe accepfance of fhe 
ALICE defecfor. Lig. ^ (lefl hand side) presenfs predictions for fhe ellipfic flow V 2 versus fransverse 
momenfum pj for fhe inclusive photons produced in 0-40% cenfral Pb-i-Pb collisions af y/s^N = 
2.76 TeV (solid red line) wifh fhe blue error band reflecfing fhe finite sfafisfics and fhe fheorefical 
uncertainly in fhe modeling of fhe cross seclions for fhe individual channels. Linally, fhe elliptic 
flow V 2 {pt) of direct photons produced in 0-40% cenfral Pb-i-Pb collisions af y/s^st = 2.76 TeV -as 
predicfed by fhe PHSD (solid red line) - is shown in Pig. Q (righf hand side); fhe blue error band is 
dominaled by fhe uncerfainfy in fhe modeling of fhe cross seclions for fhe individual channels. The 
lines presented in Pig. ^ will have lo be compared to fhe fulure experimenlal dafa in order to shed 
furlher lighl on fhe direcf photon sources and production mechanisms. 


5. Summary and outlook 

In fhis conlribulion we have calculated fhe momenfum speclra and fhe elliptic flow V 2 of di¬ 
recf phofons produced in Au-i-Au collisions af y/sji^ = 200 GeV and in Pb + Pb collisions af 
y/sj^ = 2.76 TeV using fhe microscopic PHSD fransporl approach. Lor phofon producfion we 
have incorporaled fhe inleracfions of quarks and gluons in fhe sfrongly inleracfing quark-gluon 
plasma (sQGP) {q + q ^ g + y and q{q)+ g ^ q{q) + 7 ), the photon production in the hadronic 
decays (tt —^ 7-I- 7 , T] —)• y+ 7 , m —)• tt -|- 7 , q' —^ p -|- 7 , 0 —)• q -|- 7 , tt -|- 7 ) as well as the inter¬ 

actions (71-1-71—)-p-|- 7 , p-l- 7 r—)- 7 r-|- 7 , and the bremsstrahlung radiation m+m/B^m + m/B + y) 
of mesons and baryons produced throughout the evolution of the collision. 

We find that the PHSD calculations reproduce the transverse momentum spectrum of direct 
photons as measured by the PHENIX Collaboration in Refs. [[I^, ’K]. The calculations reveal the 
channel decomposition of the observed direct photon spectrum and show that the photons produced 
in the QGP constitute at most about 50% of the direct photons with the rest being distributed 
among the other channels: mesonic interactions, decays of massive hadronic resonances and the 
initial hard scatterings. Our calculations demonstrate that the photon production in the QGP is 
dominated by the early phase (similar to hydrodynamic models) and is localized in the center of 
the fireball, where the collective flow is still rather low, i.e. on the 2-3 % level, only. Thus, the 
strong V 2 of direct photons - which is comparable to the hadronic V 2 - in PHSD is attributed to 
hadronic channels, i.e. to meson and baryon binary reactions. On the other hand, the strong V 2 
of the ’parent’ hadrons, in turn, stems from the interactions in the QGP via collisions and the 
partonic mean-filed potentials. Accordingly, the presence of the QGP shows up ’indirectly’ in the 
direct photon elliptic flow. In the future, we plan to (i) go beyond the soft photon approximation 
(SPA) in the calculation of the bremsstrahlung processes meson + meson —)• meson + meson + 7, 
meson + baryon —)• meson + baryon + 7 and (ii) will quantify the suppression at low pj due the 
Landau-Migdal-Pomeranchuk (LMP) effect. 
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The centrality dependence of the direct photon production the potential to further clarify the 
direct photon production mechanisms. We find a good agreement between the PHENIX measure¬ 
ments and PHSD calculations. In particular, the integrated thermal photon yield in PHSD was 
predicted to scale as with the exponent a = 1.5, which is in good agreement with the most 
recent measurement of a = 1.48 ± 0.08 ± 0.04 by the PHENIX Collaboration [^. This obser¬ 
vation supports the conclusion that the low transverse momentum direct photons have a strong 
contribution from the binary hadronic photon production sources, such as the meson -|- meson and 
meson + baryon bremsstrahlung. It will be important to investigate experimentally the scaling of 
the direct photon yield and flow with the number of participating nucleons Npart at EHC, too. 
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